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摘 要 :图 论 在 解决 运筹 学 .网 络 理论 .控制 论 等 领域 的 问题 中 显示 出 优越 性 。 特 别 是 最 短路 径 
问题 被 广泛 的 应 用 在 工程 .运输 等 方面 ,尤其 在 运筹 学 模型 中 最 短路 径 问题 已 是 不 可 缺少 的 一 种 
方法 。 因 此 就 其 运算 解法 进行 了 编译 ,并 用 VC + + 、Malab 两 种 程序 完成 其 算法 以 寻求 较 快 捷 的 


解法 。 
关键 词 :最 短路 径 算法 ;图 论 ; 运 筹 学 
中 图 分 类 号 :0157.6 


本 文 对 图 论 中 的 最 短路 径 问题 进行 分 析 , 就 
其 运算 解法 进行 了 编译 ,以 寻求 较 快 捷 的 解法 。 
而 关于 它 的 解法 自然 也 受到 同行 业 的 注视 , 现 就 
其 解法 从 两 种 情况 进行 介绍 : 


1 某 项 点 到 其 余 顶 点 最 短路 径 问 题 


给 定 一 个 有 向 图 C = (了 ,已 ,有 ) 以 及 某 个 项 
点 voEV( C6), 求 到 其 他 各 顶点 的 最 短路 径 的 算 
法 叫做 Dijkstra 算法 。 

在 研究 从 某 源 点 到 其 余 顶 点 最 短路 径 的 问题 时 
有 以 下 性 质 :(1) 假 设 $ 是 唯一 求 得 最 短路 径 的 终点 
的 集合 (5 的 初 态 为 空 集 ), 则 下 一 条 长 度 较 长 的 最 
短路 径 ( 设 它 的 终点 为 x) 或 者 是 弧 (v ,x); 或 者 是 中 
间 只 经 过 5 集合 中 的 顶点 ,最 后 到 达 顶 点 x 的 路 径 。 

利用 此 性 质 (2) , 按 路 径 长 度 不 减 的 顺序 产生 最 
终 路 径 。 为 此 ,引进 一 个 辅助 向 量 dist。 对 于 不 在 $ 
中 的 顶点 w, 令 dist(w) 表 示 从 wm 开始 , 且 通 过 5 中 
的 顶点 到 达 每 个 终点 w 的 最 短路 径 的 长 度 。 若 从 
vo 到 w 有 路 径 , 则 dist(w) 为 边 的 权 值 ;车 从 vo 到 w 
没有 路 径 , 则 dist(w) 为 w。 于 是 所 产生 的 下 一 条 最 
短路 径 的 终点 ,必定 是 从 w 到 所 有 的 那些 不 在 $ 
中 , 且 路 径 的 长 度 为 最 小 的 顶点 。 也 就 是 说 ,如 果 假 
设 下 一 条 最 短路 径 的 终点 是 4, 必 有 

dist(u) = minldist(w)lw¢ S$,wEV(GC)! 

由 于 确定 了 w - w 下 一 条 最 短路 径 , 于 是 把 
顶点 w 加 进 集合 $, 从 wm 开始 到 (其 中 间 点 在 5 
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中 ,而 终点 w 不 在 $ 中 ) 的 那些 路 径 , 因 为 选择 了 
u,v 便 在 S$ 中 了 ,于 是 这 条 路 径 可 以 从 wm 经 过 zx 
到 ww, 而 且 完 全 有 可 能 从 wm 经 过 v 到达 w, 比 从 
ww 不 经 过 4 到 达 w 的 路 径 要 短 , 这 时 要 修改 dist 
(w), 即 若 
dist(w) + cost(u,w) < dist(w) 
则 修改 dist(w) 为 
dist(w) = dist(wu) + cost(u,w) 
注 :cost 表示 图 的 邻接 矩阵 ,对 有 向 图 的 顶点 
从 1 至 ”编号 。 
根据 以 上 分 析 最 短路 径 算法 ,现在 分 别 介绍 
用 VC + +、Matlab 编写 程序 实现 其 算法 : 
(1) 在 VC + + 环境 中 编译 的 程序 如 下 吕 : 
(有 权 无 向 图 ) 
# include (iostream.h» 
# define MAX 50 
# define up 50000 
int cost[ MAX |] [MAX]; int dist[MAX] ,ni 
struct 
| 
int num; int pnode[ MAX |]; 
| path[ MAX]; 
void creatgraph( ) 
| 
int i, j, s, p, 1l, contin=1; 
cout《《" 顶 点 个 数 ;"; 
for (i=0;i<n;it+ +) 
| 
for (j=0;j<n;j+ +) 
cost[i][j] = cost[j][i] = up; cost[ij[i] =0; 


cin Dn; 


void shortdjs( ) 
| 
int sLMAX]; int mindis, dis, i, j, YO=0, u; 
for (i+0;i<n;i+ +) 
: 
dist[i] = cost[ VO [i}]; path[i].pnode[0] = Yo0; 
path[i] .num = 0; 
s[i] =0; 
| 
s[V0O]=1; 


for (i=1l;i<n;yi+ +) 


eee 


void dispath( ) 
| 

int 1, J; 

cout 《4”"\ n 从 V0 到 各 顶点 的 最 短路 径 长 度 如 
下 : \n’”; 

cout《《“\ t( 起 点 ~ > 终点 ) 
路 径 \ nn; 

cout (人 \t------ -=-- 


mn 
-- -------- An ; 


最 短 长 度 最 短 


ee 


void main( ) 

: 
creatgraph( ) ; shortdjs( ) ; dispath( ) ; 

| 

% 例 1 


用 VC+ + 程序 执行 结果 如 下 :; 
从 V。 到 各 顶点 的 最 短路 径 长 度 如 下 : 
起 点 一 终点 ”最 短 长 度 最短 路径 


Vo>V 2 Vo 
WT> 1 Vo ,也 
Wb 3 Vo Vi, bs 
Vo>V, 6 Vo, VV 
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(2) 在 Matlab 环境 中 编译 程序 如 (3):( 有 权 有 
向 图 ) 
建立 minroad.M 函数 
function [P,SCR)1] = minroad(i,m,W) 
% ;为 最 短路 径 的 起 始点 ,m 为 图 顶点 数 ,多 为 
图 的 带 权 邻接 惩 阵 ， 
% 不 构成 边 的 两 顶点 之 间 的 权 用 inf 表示 。 
% 3 是 矩阵 ,$ 的 每 一 列 从 上 到 下 记录 了 从 始 
点 到 终点 的 最 短路 径 所 经 顶点 的 序号 ; 
% RR 是 一 行 向 量 ,记录 了 5S 中 所 示 路 径 的 大 
小 ; 
% ”PP 是 矩阵 ,P 中 的 每 一 列 从 上 到 倒数 第 二 行 
记录 了 从 始点 到 终点 的 最 短路 径 所 经 顶点 的 序 
号 ;最 后 一 行 记录 了 该 列 所 示 路 径 的 大 小 ; 
% ”SS( 有 R) 的 意思 是 :该 位 置 根据 需要 书写 5 或 R。 
d=[];pp=[]; n=[]; n(1,1)=i; VY=1l:m; V 
(i) =[]; d= [0;i]; 
% d 的 第 二 行 是 每 次 求 出 的 最 短路 径 的 终点 ， 
第 一 行 是 最 短路 径 的 值 
p=2;la,c| = size(d); 
while ~ isempty(V) 
[td,j] = min(W(i, VY));t = Y(j); 
for k=2:c 
fa ,ji = min(d(1,k) + W(d(2,k), V)); 
2=V(jj);pp(k -1,:)= [t,t2,j1; 
end 
t= [td,t,j;ppl; [i3,14] = min(t(; ,1)); 
13= =td,n(1:2,p) = [iit(t4,2)]; 
else t5 = find(n(:,t4) ~ =0);t6 = length(t5); 
if d(2,14) = =n(t6,t4) 
n(1:t:6 +1,p) = [n(t5,t4);t(t4,2)]; 
else,n(1:3,p) = [i;d(2,t4);t(t4,2)]1; 
end;end 
d= [d,[t3;t(t4,2)]];V(t(t4,3)) =[]; 
[a,c] = size(d);p=p+1; 
end; S=n; R=4d(]1,:); 
P=[S:R); 
% 例 2 


> w=inf* ones(6); w(1,3)= 15; w(1,5) = 
20; w(1,6) =90; w(2,3)=5; w(3,4) = 50; w 
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(4,6) = 10; 

> w(5,4) =20; w(5,6)=50; w(1,2) = 1; w 
(2,5) =5; i=1; P= minroad(i,6,w) 

执行 结果 :PP = 


OOoor- 
OoOOoOLrF- 
OO 
QOD 
人 一 上 ID 一 
人 一 


0 
2 每 一 对 顶点 间 的 最 短路 径 问题 


每 次 以 一 个 顶点 为 起 点 ,重复 执行 算法 n 
次 ,可 求 得 每 一 对 顶点 之 间 的 最 短路 径 。 

求 最 短路 径 的 Floyed 算法 ,相对 较 简单 。 从 
有 向 图 的 带 树 的 邻接 矩阵 cost 出 发 ,对 有 向 图 的 
n 个 顶点 编号 。 若 从 到 j 有 弧 (i=1,2,A,n,j 
=1,2,A,n), 则 从 二 到 7 存在 一 条 长 度 为 cost (i， 
站 的 路 线 。 但 该 路 径 不 一 定 是 最 短路 径 , 尚 需 修 
改 , 修 改 的 方法 是 进行 n 次 试探 。 首 先 ,考虑 路 
和 经 (i,1,7)( 即 在 i,j 中 插 进 点 1) ,看 《i,1),《1,j) 
是 否 存 在 , 若 存 在 ,再 比较 (i,j) 与 (i,1,j) 这 两 条 
路 经 , 取 长 度 较 短 者 为 当前 求 得 的 最 短路 径 。 于 
是 这 条 求 得 的 最 短路 径 的 中 间 点 序号 不 大 于 1。 
再 在 各 对 点 i,j 中 插 进 一 个 点 2, 看 < i,…,2>， 
<2,…,j> 是 否 存在 , 若 不 存在 ,那么 当前 的 最 短 
路 径 仍 是 上 次 求 得 的 中 间 点 序号 不 大 于 1 的 最 短 
路 径 ; 若 存在 , 则 将 (i, A ,2,A,j) 的 路 径 与 前 次 求 
得 的 中 间 点 序号 不 大 于 1 的 最 短路 径 进行 比较 ， 
取 长 度 较 短 者 为 当前 最 短路 径 。 这 样 ,这 次 求 得 
的 最 短路 径 的 中 间 点 序号 不 大 于 2…… 依次 类 
推 ,直至 求 得 从 i 到 j 的 最 短路 径 。 

为 采用 邻接 矩阵 来 表示 有 向 图 , 且 每 次 修改 都 
是 在 n 阶 矩 阵 上 进行 ,首先 定义 一 个 矩阵 序列 ; 


[Se 
CN 
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A®™ ,A ,A,A'™® A 
其 中 ,4'”= cost(i,j) ,4 是 从 二 到 j 的 中 间 点 
序号 不 大 于 上 的 最 短路 径 的 长 度 , 4 的 递 推 公 
式 为 :4 = min|(AD(2,7)) ,AD ,i,k)+ 
A-DO(EII (Eken) 
{算法 分 析 】 最 短路 径 (cosi ,S(i,j),R(i,j)) 
cost 为 有 向 图 的 邻接 矩阵 ,有 向 图 的 顶点 从 
1 到 n 编号 ,S 和 RR 均 为 n 阶 方 阵 ,S( 门 为 从 - 
到 j 的 最 短路 径 长 度 , R(i, 门 为 相应 的 路 径 
循环 i 以 1 为 步 长 从 1 至 执行 
[S(i,j)<cost (i,)) 
若 izj 和 且 SS(i,j) < max 
则 RGi, 门 夺 () + ()) 
循环 以 1 为 步 长 从 1 至 nn 执行 
循环 守 以 1 为 步 长 从 1 至 n 执行 
循环 j 以 1 为 步 长 从 1 至 mn 执行 
车 S(i,j)+ S(i,j) < S(i,)) 
则 [S(i,) SC(i,k) + S(k,)) 
R(i,)))R(i,k)+ R(E,))] 
{算法 结束 1 
该 程序 的 编译 略 。 
上 述 对 最 短路 径 问 题 的 运算 程序 分 析 及 编 
译 ,在 大 量 案例 运算 中 相对 于 同类 其 他 编程 ,执行 
速度 比较 快捷 ,准确 。 该 程序 是 值得 推广 的 。 
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Solution of the shortcut problem in graph algorithms 


XIANG Rong-wu LIU Yan- jie 
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Abstract: Graph theory was a new subject which was rising quickly and was applied widely in recent years. It had 


advantages in solving the problems in many fields such as operational research, net theory and controlling theory. 


This article analyzes the shortest path problem of graph theory which come forth in engineering and transportation, 


especially about operational research model. So that by algorithm the shortest path resolution software was compiled . 


Its algorithm was completed with the two Languages, VC + + and Matlab, to seek the quicker algorithm. 


Keywords: shortest path algorithm ;graph theory ; operational research 


